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, — i ABSTRACT 

j ^ )Ve report calculations of energy levels, radiative rates, and electron impact excitation rates for transitions in Li-like ions with 12 
< Z < 20. The GRASP (general-purpose relativistic atomic structure package) is adopted for calculating energy levels and 
radiative rates, while for determining the collision strengths and subsequently the excitation rates, the Dirac atomic R-matrix 

i pode (DARC) is used. Oscillator strengths, radiative rates, and line strengths are reported for all El, E2, Ml, and M2 transitions 

^ among the lowest 24 levels of the Li-like ions considered. Collision strengths have been averaged over a Maxwellian velocity 
1ZZ4 distribution, and the effective collision strengths obtained are reported over a wide temperature range up to 10 74 K. 
^ Additionally, lifetimes are also listed for all calculated levels of the above ions. Finally, extensive comparisons are made with 
p| available results in the literature, as well as with our parallel calculations for all parameters with the Flexible Atomic Code 
+3 [fac) in order to assess the accuracy of the reported results. 
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1. Introduction 



Emission lines of Li-like ions have been observed in the spectra of a variety of astrophysical |[TJ and laboratory plasmas 
l2l-|l4). Similarly, high gain lengths have been measured in lasing plasmas, particularly for the 3d^f transitions of Al XI and 
Si XII [5j in the soft x-ray region. More recently, amplified spontaneous emission signals of 3d^-d and 3d^4f transitions of Al 
XI have been observed by Yamaguchi et al. [6|. In fact, 3d^4f (apart from other 2-3 and 3-4) transitions of Li-like ions are of 
very high interest for x-ray lasers, as demonstrated by Brown et al. [7| for ions of P, S, CI, and K. Transitions of many Li-like 
ions are also useful for temperature and density diagnostics of solar plasmas J2j, GO as well as for the determination of chemical 
abundances [9|. However to analyse observations, atomic data are required for a variety of parameters, such as energy levels, 
radiative rates (A-values), and excitation rates or equivalently the effective collision strengths (T), which are obtained from the 
electron impact collision strengths (£1). These data are also required for the modelling of fusion plasmas iflOl as many elements 
are present as impurities in the reactor walls. With this in view we have already reported calculations for Li-like ions with Z < 
1 1 ifTTl -lfLZI and in this paper we report similar results for ions with 12 < Z < 20. It may also be noted that data for transitions 
of Si XII have already been reported [13] and are therefore not included in this paper. 

Experimentally, energy levels for Li-like ions have been compiled by NIST (National Institute of Standards and Tech- 
nology) and are available at their website [http : / /www . nist . gov/pml/data/ asd . cfm| Theoretically, energy levels and 
A-values have been calculated by Nahar [14] for transitions of 15 Li-like ions in the 6 < Z < 28 range. However, collisional 
atomic data are limited to transitions among the lowest three levels of Li-like ions with 8 < Z < 92, and from these levels to 
the higher excited levels of the n < 5 configurations [15|. The calculations of Zhang et al. 03 1 are based on the relativistic 
distorted-wave (DW) method, and collision strengths (f2) are reported only at six energies above thresholds. Since resonances in 
the thresholds region have not been resolved, the subsequent results for effective collision strengths (T) are significantly under- 
estimated over a wide temperature range for a large number of transitions (in particular forbidden ones), as already demonstrated 
in our earlier papers for ions with Z < 1 1 |[TT1l - lfT2l and Si XII fT3l . 

In this paper we report a complete set of results (namely energy levels, radiative rates, and effective collision strengths) 
for all transitions among the lowest 24 levels of Li-like ions with 12 < Z < 20. Additionally, we also provide the A-values for 
four types of transitions, namely electric dipole (El), electric quadrupole (E2), magnetic dipole (Ml), and magnetic quadrupole 
(M2), because these are also required for plasma modelling. For our calculations we employ the fully relativistic GRASP (general- 
purpose relativistic atomic structure package) code for the determination of wavefunctions, originally developed by Grant et al. 
Ifl6ll and revised by Dr. P.H. Norrington. It is a fully relativistic code, and is based on the jj coupling scheme. Further rel- 
ativistic corrections arising from the Breit interaction and QED effects (vacuum polarization and Lamb shift) have also been 
included. In calculating energy levels, we have used the option of extended average level (EAL), in which a weighted (propor- 
tional to 2j+l) trace of the Hamiltonian matrix is minimized. This produces a compromise set of orbitals describing closely 
lying states with moderate accuracy. Similarly, for our calculations of fi, we have adopted the Dirac atomic R-matrix code 
(DARC) of P.H. Norrington and LP Grant (private communication). Finally, to make comparisons and to assess the accuracy of 
our results, we have performed parallel calculations from the Flexible Atomic Code (FAC) of Gu ifTTIl . available from the web- 
site http : / / sprg . ssl . berkeley . edu/ ~mf gu/ f ac/. This is also a fully relativistic code which provides a variety of 
atomic parameters, and yields results comparable to GRASP and DARC - see, for example, Aggarwal et al. [ 18] and references 
therein. Generally, there is no major discrepancy in the determination of energy levels and radiative rates from GRASP and FAC, 
but the former code allows more flexibility in terms of optimisations. Furthermore, the identification of levels belonging to a 
particular configuration is a greater problem in FAC than in GRASP, although it should not apply to simple Li-like ions. On the 
other hand, differences in collision strengths can sometimes be larger, particularly for the forbidden transitions, as discussed and 
demonstrated in some of our earlier papers (TT)-[13'|. Similarly, results for T can be significantly underestimated by the FAC 
code, particularly for forbidden transitions and at lower temperatures, due to the neglect of the contribution of resonances. Nev- 
ertheless, results from FAC will be helpful in assessing the accuracy of our energy levels, radiative rates, and collision strengths, 
and in estimating the contribution of resonances for the determination of T values. 



2. Energy levels 



The Is 2 n£ (2 < n < 5) configurations of Li-like ions give rise to the lowest 24 energy levels listed in Tables 1 (a-h). 
Included in these tables are our level energies from GRASP, obtained without and with the inclusion of Breit and QED effects, 
plus the experimental energies compiled by NIST. Also included in these tables are the theoretical energies obtained from our 
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parallel calculations from FAC and those by Nahar [14] from the Breit-Pauli R-matrix (BPRM) code of Berrington et al. Ifl9l . 
The inclusion of the Breit and QED corrections decreases the level energies (except for level 2, i.e. Is 2 2 PJ/ 2 ) by a maximum of 
~ 0.03 Ryd, depending on the ion. As the ion nuclear charge Z increases, so does the effect of the Breit and QED corrections. 
For example, energy levels of Mg X decrease by ~ 0.005 Ryd, but for Ca XVIII by ~ 0.03 Ryd. Clearly, the importance of 
relativistic effects increases with increasing Z. Our Breit and QED corrected energies from GRASP agree closely (generally 
within 0.01 Ryd) with those of NIST and the orderings are also nearly the same, although there are some noticeable differences, 
such as for the 4d,5d 2 D 3 / 25 / 2 levels of K XVII and Ca XVIII. Degeneracy for these levels in the NIST compilations is up to 
0.061 Ryd, whereas it is ~ 0.01 Ryd in all the calculations listed in Tables 1 (a-h). Furthermore, the NIST energies for some of 
the levels are non-degenerate - see, for example, the 4p,5p 2 P°/ 2 3/2 l eve l s of Mg X and 5d 2 D 3 / 2 ,5/2 levels of P XIII. Similarly, 
NIST energies are not available for some of the levels - see, for example, the 5f 2 F^ 2 7 , 2 and 5g 2 G- 7 / 2 9 / 2 levels of P XIII, CI 
XV, Ar XVI, K XVII, and Ca XVIII. 

The energies obtained from FAC are generally comparable with the experimental compilations of NIST and/or our cal- 
culations from GRASP. However, the orderings of some of the levels are slightly different from FAC than from the experimental 
and other theoretical energies. As examples, see 12-13 and 19-24 levels of Mg X and 16-24 of Ca XVIII. Nevertheless, energy 
differences among these levels are very small. Finally, the BPRM energies of Nahar [14| have the same orderings as NIST and 
our calculations from GRASP. The BPRM energies also generally agree in magnitude, except for 4s,5s 2 Si/ 2 levels for which 
they are lower by up to 0.2 Ryd, depending on the ion. Overall, we may state that there is no (major) discrepancy between theory 
and experiment for the energy levels of Li-like ions. 



3. Radiative rates 



The absorption oscillator strength (Jij) and radiative rate Aji (in s 1 ) for a transition i — > j are related by the following 
expression: 



fa = — r^ A "<-- 1 - = L49 x io^a^-M)^ (i) 

07T e LOi 

where m and e are the electron mass and charge, respectively, c is the velocity of light, Xji is the transition energy/wavelength in 
A, and LUi and uij are the statistical weights of the lower (i) and upper (j) levels, respectively. Similarly, the oscillator strength 
(dimensionless) and the line strength S (in atomic units, 1 a.u. = 6.460 x 10 -36 cm 2 esu 2 ) are related by the following standard 
equations: 



for the electric dipole (El) transitions 



for the magnetic dipole (Ml) transitions 



A,, = 2 - 6974 * 10 'V and /y = i21^V, ,3, 



for the electric quadrupole (E2) transitions 



UJ 3 A ji A yi^i 
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and for the magnetic quadrupole (M2) transitions 



1.4910 > ;M » SM2 and (5) 

In Tables 2 (a-h) we present transition energies/wavelengths (A, in A), radiative rates (A^, in s -1 ), oscillator strengths 
(fij, dimensionless), and line strengths (S, in a.u.), in length form only, for all 92 electric dipole (El) transitions among the 24 
levels of Li-like ions with 12 < Z < 20. The indices used to represent the lower and upper levels of a transition have already 
been defined in Tables 1 (a-h). Similarly, there are 100 electric quadrupole (E2), 79 magnetic dipole (Ml), and 108 magnetic 
quadrupole (M2) transitions among the 24 levels. However, for these transitions only the A-values are listed in Tables 2 (a-h), 
and the corresponding results for f- or S-values may be easily obtained using Eqs. (1-5). 

In Tables A-H we compare our A-values from the GRASP and FAC codes for some of the common El transitions with 
those of Nahar [14| from the BPRM code of Berrington et al. [ 19 1. Also included in these tables are the f-values from our GRASP 
calculations. In general, the A-values from the GRASP and FAC codes agree satisfactorily for a majority of transitions, although 
there are differences of up to 30% for a few, such as 1-17 and 1-18 of Mg X, for which there are no discrepancies with the 
BPRM calculations of Nahar. However, these differences between two independent calculations decrease with increasing Z, i.e. 
reduce to 13% for the corresponding transitions of Ca XVIII. On the other hand, there are significant differences with the Nahar 
calculations, of up to two orders of magnitude, particularly for the comparatively weaker (f < 0.01) transitions, such as 7-11, 
7-18, and 12-18. For these transitions the BPRM A-values of Nahar are lower for all Li-like ions, as also noted for ions with Z < 
11 [ 1 2 1 - [ 1 3 1 . The A-values for such weak transitions from the BPRM code are likely to be less accurate, as discussed by Hibbert 
||20l and Aggarwal et al. ED . In addition, several transitions are missing from the Nahar calculations, and examples include: 
2p 2 P° /2 - 3s 2 S 1/2 (2-4), 2p 2 P° /2 - 3s 2 S 1/2 (3^4), and 2p 2 P° /2 - 3d 2 D 3/2 (3-7). Finally, there are differences of up to a 
factor of four between the BPRM A-values of Nahar and our calculations from the GRASP and FAC codes for several transitions, 
particularly those among levels 16 and higher. For all such transitions, the energy differences (A-Ey) are very small and hence 
small change in these affect the A-values significantly. 

One of the general criteria to assess the accuracy of radiative rates is to compare the length and velocity forms of the 
f- or A-values. However, such comparisons are only desirable, but are not a fully sufficient test to assess accuracy, as different 
calculations (or combinations of configurations) may give comparable f-values in the two forms, but entirely different results in 
magnitude. Generally, there is good agreement between the length and velocity forms of the f-values for strong transitions, but 
differences between the two forms can sometimes be substantial even for some very strong transitions, as demonstrated through 
various examples by Aggarwal et al. [22 J . However, for transitions of Li-like ions, there is no discrepancy between the two forms, 
and this is irrespective of the magnitude of a transition. Furthermore, to assess the importance of configuration interaction (CI), 
we have performed another calculation including the levels of the n = 6 configurations. Only for two transitions, namely 2-16 
(2p 2 P° /2 - 5s 2 Si /2 ) and 3-16 (2p 2 P° /2 - 5s 2 Si /2 ), have the A-values decreased by 30% for Mg X and by 20% for Ca XVIII. 
For the remaining transitions the effect of additional CI is much less noticeable. Based on this and other comparisons discussed 
above, we may state that the A-values listed in Tables 2 (a-h) from our GRASP calculations are accurate to better than 20% for 
all transitions and all ions, and hence should be the best currently available. 



4. Lifetimes 



The lifetime r for a level j is defined as follows: 



Ei A 3 



(6) 



Since this is a measurable parameter, it provides a check on the accuracy of the calculations. Therefore, in Tables 1 
(a-h) we have also listed our calculated lifetimes from the GRASP code, which include the contributions from four types of 
transitions, i.e. El, E2, Ml, and M2. To our knowledge, no measurements are available for lifetimes in Li-like ions with 12 < Z 
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< 20. However, Ishii et al. |23l have measured lifetimes from beam foil experiments for the 2p 2 P^ 2 3/2 l eve l s °f CI XV, namely 
1 .0±0. 1 and 0.71 ±0.03 ns, respectively. Both these measured values agree very well with our calculated lifetimes. Froese-Fischer 
has calculated lifetimes for all the Li-like ions under discussion here and has posted her results for some levels up to 4s 2 Si/ 2 at 
her website: http : / /www . vuse . vanderbilt . edu/ ~ cff/mchf_col lection/. There are no discrepancies for any 
of the levels between our and her calculations. We hope the present results for other levels will be useful for future comparisons 
and may also encourage experimentalists to measure lifetimes for the Li-like ions. 



5. Collision strengths 



For the computation of collision strengths (ST) we have employed the Dirac atomic R-matrix code (DARC), which includes 
the relativistic effects in a systematic way, in both the target description and the scattering model. It is based on the jj coupling 
scheme, and uses the Dirac-Coulomb Hamiltonian in the i?-matrix approach. The i?-matrix radii adopted for Li-like ions with 12 
< Z < 20 (except Z = 14) are 9.60, 8.64, 7.20, 6.40, 5.92, 5.44, 5.28, and 5.12 au, respectively. For all eight ions, 56 continuum 
orbitals have been included for each channel angular momentum for the expansion of the wavefunction. This allows us to compute 
SI up to an energy of 180, 215, 310, 390, 450, 540, 570, and 600 Ryd, for 12 < Z < 20, respectively. Furthermore, these energy 
ranges are sufficient to calculate values of T up to T e = 10 7 4 K, i.e. well above the temperature of maximum abundance of 
these ions in ionisation equilibrium - see Bryans et al. Il24l . The maximum number of channels for a partial wave is 108, and 
the corresponding size of the Hamiltonian matrix is 6086. To obtain convergence of SI for all transitions and at all energies, we 
have included all partial waves with angular momentum J < 60, although a larger number would have been preferable for the 
convergence of some allowed transitions, especially at higher energies. However, to account for neglected higher partial waves, 
we have included a top-up, based on the Coulomb-Bethe [25 1 approximation for allowed transitions and geometric series for 
others. 

For illustration, in Figs. 1-3 we show the variation of SI with angular momentum J for three transitions of Ar XVI, namely 
1-3 (2s 2 Si /2 - 2p 2 P3 /2 ), 1-8 (2s 2 Si /2 - 3d 2 D 5/2 ), and 1-10 (2s 2 Si /2 - 4p 2 V° 1/2 ), respectively, and at five energies of 100, 
200, 300, 400, and 500 Ryd. Values of SI have fully converged for all forbidden (such as 1-8) and most of the allowed (such as 
1-10) transitions. However, for some allowed transitions, such as 1-3 shown in Fig. 1, our adopted range of partial waves (J < 
60) is not sufficient for the convergence of il at most of the energies, particularly above thresholds. For such transitions a top-up 
has been included as mentioned above, and has been found to be appreciable. A similar conclusion applies to the transitions in 
other Li-like ions. 

In Tables 3 (a-h) we list our values of SI for resonance transitions in Li-like ions with 12 < Z < 20 at energies above 
thresholds. The indices used to represent the levels of a transition have already been defined in Tables 1 (a-h). Very limited 
data are available for comparison, as already stated in section 1 . Therefore, to make an accuracy assessment of the values of SI, 
we have performed another calculation using the FAC code of Gu ifTTl . This code is also fully relativistic, and is based on the 
well-known and widely-used distorted-wave (DW) method. Furthermore, the same CI is included in FAC as in the calculations 
from DARC. Therefore, also included in Tables 3 (a-h) for comparison purposes are the SI values from FAC at a single excited 
energy E jt which corresponds to ~ 140, 170, 230, 400, 450, 510, 580, and 450 Ryd for 12 < Z < 20, respectively. 

Generally the two sets of SI from DARC and FAC agree to within 20% for a majority of transitions. However for ~ 15% 
of the transitions there are differences of over 20%, and for ^5% the discrepancies between the two independent calculations 
are over 50%. Most of the transitions for which the discrepancies are large are forbidden, and one of the main reasons is the 
anomalous behaviour of SI in the calculations from the FAC code. To illustrate this, in Fig. 4 we compare the two sets of S7 from 
DARC and FAC for three allowed transitions of Ar XVI, namely 4-6 (3s 2 S 1/2 - 3p 2 ?3 /2 ), 6-8 (3p 2 P° /2 - 3d 2 D 5/2 ), and 10-12 
(4p 2 P°/ 2 - 4d 2 D 3 / 2 ). For these three and many other allowed transitions there are no major discrepancies between the two sets 
of SI, although there is a fall in magnitude for the 10-12 (and some other) transitions at the higher end of the energy range for all 
ions. In Fig. 5 we show similar comparisons for three forbidden transitions, namely 6-1 1 (3p 2 Py 2 ~ 2 ^3/2^ 8—13 (3d 2 D 5 / 2 
-4d 2 D 5 / 2 ), and 10-13 (4p 2 P^ 2 -4d 2 D 5 / 2 ). Differences for these three (and some other) transitions are generally within 50%, 
but are up to a factor of two for a few, such as 12-16, 13-16, 16-18, 19-21, and 20-22. For some transitions, such as 6-11, there 
are differences over the entire energy range, while for others such as 8-13, the discrepancies are mainly at the lower end of the 
energy range, and completely disappear at large energies. In some cases, such as 10-13, the differences between the two sets of SI 
are mainly at the higher end of the energy range due to a sudden anomaly in S7 from the FAC code. Discrepancies in the values of ft 
for some transitions from DARC and FAC are partly due to the differences in the methodology, i.e. i?-matrix and DW, respectively, 
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and partly because FAC is designed to generate a large amount of atomic data in a comparatively very short period of time, and 
without too much loss of accuracy. This aim is achieved by employing various interpolation and extrapolation techniques, which 
result in significant differences in the values of ft for some transitions and at some energies, as noted above. Overall, we may 
state that for transitions for which there is no anomaly in the calculations of 57 from the FAC code, the differences between the 
DARC and FAC calculations of collision strengths are generally within 20%. However, the discrepancies are significant for those 
transitions for which values of £1 from the FAC code are anomalous. A similar conclusion applies to transitions in other Li-like 
ions. Finally, based on the comparisons made above as well as for other Li-like ions |[TTI - |[T3l . we estimate our results for ft to 
be accurate to better than 20% for a majority of the transitions of all eight Li-like ions. 



6. Effective collision strengths 



Excitation rates, along with energy levels and radiative rates, are required for plasma modelling, and are determined from 
the collision strengths (SI). Since the threshold energy region is dominated by numerous closed-channel (Feshbach) resonances, 
values of ft need to be calculated in a fine energy mesh in order to accurately account for their contribution. Furthermore, in a 
hot plasma electrons have a wide distribution of velocities, and therefore values of f2 are generally averaged over a Maxwellian 
distribution as follows: 



T(T e ) = / n(E)exp(-Ej/kT e )d(Ej/kT e ), (7) 
Jo 

where k is Boltzmann constant, T e is electron temperature in K, and Ej is the electron energy with respect to the final (excited) 
state. Once the value of T is known the corresponding results for the excitation q(i, j) and de-excitation q(j, i) rates can be easily 
obtained from the following equations: 



8.63 x 10- 6 , 



q(i,j) = j—Texpt-Ea/fciy crn^s" 1 (8) 

UJiT e 



and 

8.63xl0- 6 „ .. 



?0M) = - T7^ T cms > < 9 > 

UljT e 

where and uij are the statistical weights of the initial (i) and final (j) states, respectively, and Eij is the transition energy. 
The contribution of resonances may enhance the values of T over those of the background values of collision strengths (£Ib), 
especially for the forbidden transitions, by up to a factor of ten (or even more) depending on the transition and/or the temperature. 
Similarly, values of need to be calculated over a wide energy range (above thresholds) in order to obtain convergence of the 
integral in Eq. (7), as demonstrated in Fig. 7 of Aggarwal and Keenan [26]. 



To delineate resonances, we have performed our calculations of il in the threshold region at over 35000 energies, depend- 
ing on the ion. Close to thresholds (~ 0.1 Ryd above a threshold) the energy mesh is 0.001 Ryd, and away from thresholds is 
0.002 Ryd. Thus care has been taken to include as many resonances as possible, and with as fine a resolution as is computation- 
ally feasible. The density and importance of resonances can be appreciated from Figs. 6-8, where we show in the thresholds 
region for three transitions of Ar XVI, namely 1-2 (2s 2 S 1/2 - 2p 2 V° l/2 ), 1-3 (2s 2 S 1/2 - 2p 2 ?l /2 ), and 2-3 (2p 2 P° /2 - 2p 
2 P^ 2 ). Similar resonances are observed for transitions in other Li-like ions. These resonances affect the values of T particularly 
towards the lower end of the temperature range. 



Our calculated values of T are listed in Tables 4 (a-h) over a wide temperature range up to 10 7 4 K, suitable for applications 
in solar and other plasmas. As stated in section 1, there are no similar results available for comparison purposes for most of the 
ions. Therefore, as for collision strengths we have also calculated values of T from our non-resonant VL data from the FAC code, 
which are included in Tables 4 (a-h) at the lowest and the highest calculated temperatures for each ion. This provides a ready 
comparison between the two independent calculations, and also gives an estimate of the importance of resonances. For example, 
at T e = 10 5 5 K, our resonances-resolved values of T differ by over 20% for 120 (~ 40%) transitions of Mg X, and in some 
instances by up to an order of magnitude, see for example: 2-4, 3-4, and 7-9. Many of these are resonance transitions, such as 
1-5/6/7/8/10/11, and our results from DARC are invariably higher due to the inclusion of resonances. However, at T e = 10 6 9 K, 
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only 24 (~ 9%) transitions differ by over 20%, with a maximum discrepancy of a factor of two for transitions such as: 2-17, 
5-17, and 6-18. For the 12-22 (4d 2 D 3 / 2 - 5f 2 F^ 2 ) transition, values of T from the FAC code are higher by a factor of 2.7 
because of a sudden rise in the corresponding values of il, as shown for the 10-13 transition in Fig. 5. Overall we may state that 
the contribution of resonances is more dominant and appreciable at lower temperatures. Similar conclusions apply for transitions 
in other Li-like ions. 

The only Li-like ion considered here for which results for T are available in the literature is Ar XVI. Whiteford et al. 
||27l have reported values of T for all transitions among the lowest 24 levels over a wide range of electron temperature. They 
have adopted the standard i?-matrix code of Berrington et al. [19] and have included one -body relativistic operators in the 
calculations, which should be sufficient for a moderately heavy ion such as Ar XVI. Furthermore, they included a comparable 
large range of partial waves with angular momentum J < 58 and resolved resonances in the thresholds region in a narrow energy 
mesh of ~ 0.002 Ryd. Therefore, their results for T should be comparable to our calculations from the fully relativistic DARC 
code. A comparison between the two independent sets of T is indeed highly satisfactory for a majority of transitions. For only 
10 transitions are the differences over 20%, and the discrepancies are particularly large for seven, namely 7-8, 12-13, 14-15, 
19-20, 21-22, 21-24, and 23-24. All of these transitions are forbidden, there is no (major) discrepancy between the ft values 
from the DARC and FAC codes, and there is no large variation of their fl values with increasing energy, i.e. is quite stable. 
Therefore, the differences in the values of T should not be due to the corresponding differences in the values of il, but because 
of the contribution of resonances. However, all of these transitions belong to the same n and therefore their AEij are very small, 
but resonances in collision strengths are neither too dense nor too high in magnitude. In fact, 7-8 (3d 2 D 3 / 2 - 3d 2 D 5 / 2 ) is the 
only transition for which both the thresholds range and resonances are appreciable, as shown in Fig. 9. 

In Fig. 10 we compare values of T from the BPRM calculations of Whiteford et al. Il27l with our own from the DARC 
code for three transitions, namely 7-8 (3d 2 D 3/2 - 3d 2 D 5/2 ), 12-13 (4d 2 D 3/2 - 4d 2 D 5/2 ), and 14-15 (4f 2 F° /2 - 4f 2 F° /2 ). 
At the lowest common temperature of 5.12xl0 4 K, values of T from the BPRM calculations are higher, by up to a factor 
of ~30, for the above three and some other transitions, noted above. As the electron temperature increases the discrepancy 
between the two calculations decreases significantly, and towards the higher end of the temperature range the two independent 
calculations become comparable. The resonances for the 7-8 transition are neither too high in magnitude nor dense (see Fig. 
9), and hence cannot make a large contribution to the determination of T as shown by our calculations from DARC. Clearly, 
the BPRM calculations of T have overestimated the contribution of resonances for the above noted seven transitions. The main 
reason for such a large discrepancy between the two i?-matrix calculations for a few transitions, particularly those between the 
degenerate levels of a state, is the limitation of the BPRM method, as recently discussed and demonstrated by Bautista et al. ||28l . 
The jj coupling method adopted in the DARC code includes fine-structure in the definition of channel coupling and hence gives 
a better representation of resonances, particularly for transitions within the degenerate levels of a state. These near-threshold 
resonances affect the subsequent results of T, especially at lower temperatures as shown in Fig. 10. We discuss this further 
below. 

Another source of higher values of T from the BPRM calculations of Whiteford et al. [27], for some transitions and at 
lower temperatures, could be their inclusion of higher levels of the lsnln'i' configurations. These additional levels in the 226- 
270 Ryd energy range, well above the highest threshold of our calculations, i.e. ~ 57 Ryd, may give rise to some resonances 
which may affect the values of T. However, it may be true for the calculations of T at very high temperatures, such as T e = 
10 7 4 K (equivalent to ~ 163 Ryd), but resonances will not contribute to T values at lower temperatures, such as T e = 10 47 
K (equivalent to ~ 0.32 Ryd) until there are very high spurious resonances. Therefore, we are confident of our results listed 
in Tables 4 (a-h) which are assessed to be accurate to better than ^20% for a majority of transitions and over a wide range of 
temperature. However, the contribution of higher neglected resonances, from the 1 s 2 n£ (n > 6) and 1 sn£n'£' configurations, may 
affect these results towards the higher end of the temperature range, as discussed above, and particularly for transitions involving 
the levels of the n = 5 configurations. 

7. Conclusions 



In this paper we have presented results for energy levels and radiative rates for four types of transitions (El, E2, Ml, and 
M2) among the lowest 24 levels of Li-like ions with 12 < Z < 20 belonging to the n < 5 configurations. Additionally, lifetimes 
of all the levels have been reported, although measurements or other theoretical results are not available for comparison for most 
of the ions/levels. However, based on a variety of comparisons among various calculations from the GRASP and FAC codes, our 
energy levels are assessed to be accurate to better than 1 %, and the results for radiative rates, oscillator strengths, line strengths, 
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and lifetimes are assessed to be accurate to better than 20% for a majority of strong transitions (levels). Similarly, the accuracy 
of our results for collision strengths and effective collision strengths is estimated to be better than 20% for a majority of the 
transitions. This accuracy estimate is based on the comparison made between two independent calculations performed with the 
DARC and FAC codes, as well as with other results available in the literature, particularly for Ar XVI. Additionally, we have 
considered a large range of partial waves to achieve convergence of 57 at all energies, included a wide energy range to accurately 
calculate the values of T up to T e = 10 7 4 K, and resolved resonances in a fine energy mesh to account for their contributions. 
Hence we see no obvious deficiency in our calculated results. However, the present results for effective collision strengths for 
transitions involving the levels of the n = 5 configurations may (perhaps) be improved by the inclusion of the levels of the n 
= 6 configurations, as has been noted for the A-values. Finally, we would like to note that inclusion of radiation and Auger 
damping may reduce the contribution of resonances, particularly at lower temperatures, for some of the transitions. However, it 
is certainly true for transitions involving the inner-shell, as demonstrated by Whiteford et al. (27 j, but there is no evidence of 
similar (significant) effect for the transitions considered in this paper. We believe the present set of complete results for radiative 
and excitation rates for eight Li-like ions are probably the best available todate and will be highly useful for the modelling of a 
variety of plasmas. 
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Captions for Figures 



Figure 1 Partial collision strengths for the 2s 2 Si/ 2 - 2p 2 P 3 / 2 (1-3) transition of Ar XVI, at five energies of: 

100 Ryd (circles), 200 Ryd (triangles), 300 Ryd (stars), 400 Ryd (squares), and 500 Ryd (diamonds). 
Figure 2 Partial collision strengths for the 2s 2 Si/ 2 - 3d 2 D 5 / 2 (1-8) transition of Ar XVI, at five energies of: 

100 Ryd (circles), 200 Ryd (triangles), 300 Ryd (stars), 400 Ryd (squares), and 500 Ryd (diamonds). 
Figure 3 Partial collision strengths for the 2s 2 Si/ 2 - 4p 2 P°/ 2 (1-10) transition of Ar XVI, at five energies of: 

100 Ryd (circles), 200 Ryd (triangles), 300 Ryd (stars), 400 Ryd (squares), and 500 Ryd (diamonds). 
Figure 4 Comparison of collision strengths from our calculations from DARC (continuous curves) 

and FAC (broken curves) for the 4-6 (circles: 3s 2 Si/ 2 - 3p 2 P 3 / 2 ), 6-8 (triangles: 3p 2 P 3 /, - 3d 2 D 5 / 2 ), and 

10-12 (stars: 4p 2 P° /2 - 4d 2 D 3/2 ) allowed transitions of Ar XVI. 
Figure 5 Comparison of collision strengths from our calculations from DARC (continuous curves) and 

FAC (broken curves) for the 6-1 1 (triangles: 3p 2 P° /2 - 4p 2 P 3/2 ), 8-13 (stars: 3d 2 D 5/2 - 4d 2 D 5/2 ), and 

10-13 (circles: 4p 2 P° /2 - 4d 2 D 5/2 ), forbidden transitions of Ar XVI. 
Figure 6 Collision strengths for the 2s 2 S;iy 2 - 2p 2 P°y 2 (1-2) transition of Ar XVI. 
Figure 7 Collision strengths for the 2s 2 Si/ 2 - 2p 2 P° /2 (1-3) transition of Ar XVI. 
Figure 8 Collision strengths for the 2p 2 P° /2 - 2p 2 P° /2 (2-3) transition of Ar XVI. 
Figure 9 Collision strengths for the 3d 2 D 3/2 - 3d 2 D 5/2 (7-8) transition of Ar XVI. 
Figure 10 Comparison of effective collision strengths for the 7-8 (circles: 3d 2 D 3 / 2 - 3d 2 D 5 / 2 ), 

12-13 (triangles: 4d 2 D 3/2 - 4d 2 D 5/2 ), and 14-15 (stars: 4f 2 F° /2 - 4f 2 F° /2 ) 

transitions of Ar XVI. Continuous and broken curves are from the present DARC and 

earlier BPRM codes [27|, respectively. 
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Figure 1. Partial collision strengths for the 1-3 transition of Ar XVI. 



in 




J + 1 

Figure 1: Partial collision strengths for the 2s 2 S 1 / 2 - 2p 2 ^/ 2 (1-3) transition of Ar XVI, at five energies of: 100 Ryd (circles), 
200 Ryd (triangles), 300 Ryd (stars), 400 Ryd (squares), and 500 Ryd (diamonds). 
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Figure 2. Partial collision strengths for the 1-8 transition of Ar XVI. 
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Figure 2: Partial collision strengths for the 2s 2 Si/ 2 - 3d 2 D 5 / 2 (1-8) transition of ArXVI, at five energies of: 100 Ryd (circles), 
200 Ryd (triangles), 300 Ryd (stars), 400 Ryd (squares), and 500 Ryd (diamonds). 
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Figure 3. Partial collision strengths for the 1-10 transition of Ar XVI. 
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Figure 3: Partial collision strengths for the 2s 2 S!/ 2 -4p 2 PJ/ 2 (1—10) transition of ArXVI, at five energies of: 100 Ryd (circles), 
200 Ryd (triangles), 300 Ryd (stars), 400 Ryd (squares), and 500 Ryd (diamonds). 
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Figure 4. Comparison of collision strengths for the 4-6, 6-8, and 10-12 allowed transitions of Ar XVI. 
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Figure 4: Comparison of collision strengths from our calculations from DARC (continuous curves) and FAC (broken curves) for 
the 4-6 (circles: 3s 2 S 1/2 - 3p 2 P£ /2 ), 6-8 (triangles: 3p 2 P° /2 - 3d 2 D 5/2 ), and 10-12 (stars: 4p 2 P° /2 - 4d 2 D 3/2 ) allowed 
transitions of Ar XVI. 
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Figure 5. Comparison of collision strengths for the 6—11, 8—13 and 10—13 forbidden transitions of Ar XVI. 
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Figure 5: Comparison of collision strengths from our calculations from DARC (continuous curves) and FAC (broken curves) for 
the 6-11 (triangles: 3p 2 P° /2 -4p 2 P° /2 ), 8-13 (stars: 3d 2 D 5/2 -4d 2 D 5/2 ), and 10-13 (circles: 4p 2 P° /2 -4d 2 D 5/2 ), forbidden 
transitions of Ar XVI. 
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Fig. 6. Collision strengths for the 1 -2 transition of Ar XVI 
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Figure 6: Collision strengths for the 2s 2 Si/ 2 - 2p 2 P° /2 (1-2) transition of Ar XVI. 
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Fig. 7. Collision strengths for the 1 —3 transition of Ar XVI 
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Figure 7: Collision strengths for the 2s 2 S 1/2 - 2p 2 P° /2 (1-3) transition of Ar XVI. 
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Fig. 8. Collision strengths for the 2-3 transition of Ar XVI 
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Figure 8: Collision strengths for the 2p 2 P° /2 - 2p 2 P° /2 (2-3) transition of Ar XVI. 
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Fig. 9. Collision strengths for the 7-8 transition of Ar XVI 
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Figure 9: Collision strengths for the 3d 2 D 3/2 - 3d 2 D 5/2 (7-8) transition of Ar XVI. 
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Figure 10. Comparison of effective collision strengths for the 7-8, 12 — 13 and 14—15 transitions of Ar XVI. 
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Figure 10: Comparison of effective collision strengths for the 7-8 (circles: 3d 2 D 3 / 2 - 3d 2 D 5 / 2 ), 12-13 (triangles: 4d 2 D 3 / 2 
- 4d 2 D 5 / 2 ), and 14-15 (stars: 4f 2 F^ 2 - 4f 2 Fyy 2 ) transitions of Ar XVI. Continuous and broken curves are from the present 
DARC and earlier BPRM codes |27|, respectively. 
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Table A. Comparison of radiative rates (A- values, s : ) for some transitions of Mg X. (a±b = ax 10 ±6 ). 
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Table B. Comparison of radiative rates (A- values, s x ) for some transitions of Al XI. (a±b = ax 1CP ). 
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Table C. Comparison of radiative rates (A- values, s 1 ) for some transitions of P XIII. (a±b = ax I0 ±b ). 
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GRASP: Present results from the GRASP code 
FAC: Present results from the FAC code 
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Table D. Comparison of radiative rates (A- values, s 1 ) for some transitions of S XIV. (a±b = ax 10 ±fc ). 
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Table E. Comparison of radiative rates (A- values, s x ) for some transitions of CI XV. (a±b = ax 10 ±b ). 
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GRASP: Present results from the GRASP code 
FAC: Present results from the FAC code 
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Table F. Comparison of radiative rates (A- values, s 1 ) for some transitions of Ar XVI. (a±b = ax 10 ±b ). 
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Table G. Comparison of radiative rates (A- values, s x ) for some transitions of K XVII. (a±b = ax 10 ). 
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GRASP: Present results from the GRASP code 
FAC: Present results from the FAC code 
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Table H. Comparison of radiative rates (A- values, s x ) for some transitions of Ca XVIII. (a±b = ax 10 ±6 ). 
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